LINEAR STABILITY ANALYSIS FOR PERIODIC TRAVELING WAVES OF 
THE BOUSSINESQ EQUATION AND THE KGZ SYSTEM 

SEVDZHAN HAKKAEV, MILENA STANISLAVOVA, AND ATANAS STEFANOV 

Abstract. The question for linear stability of spatially periodic waves for the Boussinesq equa- 
tion (the cases p — 2,3) and the Klein-Gordon-Zakharov system is considered. For a wide class of 
fSJ solutions, we completely and explicitly characterize their linear stability (instability respectively), 

T-H when the perturbations are taken with the same period T. In particular, our results allow us to 

completely recover the linear stability results, in the limit T — )■ oo, for the whole line case. 

(N 
X) 

G^ 1. Introduction 



In this paper we will be interested in the stability of spatially periodic waves for certain models, 
which involve second derivative in time. Our interest will be mainly in two PDE - the Boussinesq 
equation and the Klein-Gordon-Zakharov system, although the methods that we develop here 
will certainly find applications in other related models. 

The Cauchy problem for the Boussinesq equation, with periodic boundary conditions, is 



< 



' ' (1) Uu + Uxxxx-Uxx + if{u))xx=0, (t, x) G R+ X [0,T], 

^ where f{u) will for the most part be f{u) = vP,p > 1. This is a model that was derived by 

^r^ Boussinesq, for p = 2, but was subsequently studied by many authors, both in the periodic 
and whole line context. We now review the current results regarding the well-posedness properties 
of the Boussinesq equation. While we have a very satisfactory theory for the local solutions, see 
below, the the global well-posedness does not hold. More precisely, even if one requires smooth 

^ compactly supported data, the solutions may develop singularities in finite time, [8]. This makes 
the stability questions, which is the main subject of this article even more relevant and interesting. 

^— I In the whole line scenario, local well-posedness was established by Bona and Sachs, in the 

^ Sobolev spaces f/'2 + (Ri) x ff2 + (Ri)j [8]. Further contributions were made by Tsutsumi and 

k> Mathashi, |27j . Linares, |21j (who also showed global existence for small data). Farah, fl3j has 
shown well-posedness in i?*(R^) x H^~'^(R}), when s > —1/4 and the space if" is defined via 
^ = {u : Ux G H'^~^ (R})} . Kishimoto and Tsugava, [20J have finally shown well-posedness for 

all s > —1/2, which is likely to be sharp. 

Regarding the case of periodic boundary conditions, Fang and Grillakis, [12] who have estab- 
hshed local well-posedness in if*(T) x H'^~'^{T), s > (when 1 < p < 3 in ([T])). This result 
was improved for p = 2to,s> —1/4 by Farah and Scialom, [T^. Oh and Stefanov have recently 
shown local well-posedness in if^(T) x H'-^iT), s > -3/8, 
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Our other main object of investigation will be the Klein-Gordon-Zakharov system, which is 
given h:^ 

( utt -Uxx + u + uv = 0, {t, x) G X or {t, x) £ x [0, T] 
\ ntt- n^x = \{\u\'^)xx 

This system describes the interaction of a Langmuir wave and an ion sound wave in plasma. More 
precisely, u is the fast scale component of the electric field, whereas n denotes the deviation of 
ion density, |281lllj. The system ^ is locally well-posed in various function spaces, (see Ginibre- 
Tsutsumi-Velo, [l^ and Ozawa-Tsutaya-Tsutsumi, [23]). In our previous paper, [19], we have 
shown that the Cauchy problem ^ is locally well-posed (both in priodic and whole line contexts) 
in H"" X X H""'^ x whenever a > 1/2. In [23j and [24J, the authors have shown that 

small data produces solutions that persist globally, whereas large solutions are generally expected 
to blow up in finite time. 

The stability of periodic traveling waves have been studied extensively in the last decade. 
The nonlinear stability of periodic waves for Koreteweg-de Vries eqtuation based on the Jacobi 
elliptic functions of cnoidal type was considered in and for modified Korteweg-de Vries and 
nonlinear Schrodinger equations of dnoidal type in|3]. In [T7], the authors considered the nonlinear 
stability of periodic waves for generalized Benjamin-Bona-Mahony equation. Recently, Arruda 
[5] considered the nonlinear stability of periodic traveling waves for Boussinesq equation. The 
approach is based on the theory developed in [6l [7] and [16] for the stability of solitary waves. 

In this paper, we consider the linearized stability of periodic traveling waves for the Boussinesq 
model (the cases p = 2,3) and the Klein-Gordon-Zakharov system. These are the cases of 
second order in time models, for which we can explicitly write the solutions in elliptic functions 
(and moreover, we can explicitly compute the relevant portion of the spectrum of the linearized 
operators). While this certainly helps in the analysis, we believe that our results should be 
generalizable to all values of p > 1. The main tool is the theory for linearized stability for such 
models, developed recently by the second and third author, [26] . 

The paper is organized as follows. In Section [2| we present the construction of our main object 
of study - the periodic traveling waves. This is not a new material by any means, but we do 
it in order to single out the solutions of interest^ and to introduce some notations. In Section 
[3] we setup the linear stability problem, after which we present the main results. In Section [4] 
we outline the theory for linearized stability for second order ODE in [5H], and point out to the 
relevant spectral theoretic results about their linearized operators. In Section [5j we prove the 
main results - Theorems [l| [2| about the Boussinesq model, while in Section [6j we prove Theorem 
[3] about the KGZ system. 



2. Construction of the periodic traveling waves 



In this section, we show a glimpse of the construction of the periodic waves - in Section 2.1 for 



the Boussinesq equation (when p = 2,3) and in Section 2.2 for the KGZ system. 



2.1. Construction of the traveling wave solutions for the Boussinesq equation. Apply- 
ing the traveling wave ansatz, one sees that there is an one-parameter family of traveling waves of 



"'^The coefficient | in front of tfie non-linear term (|up)a;a; is non-standard, but ratlier adopted for convenience of 
our presentation. In particular, it fielps create a self-adjoint linearized operator, wliicli otherwise can be achieved 
via a simple change of the time variable. 

^note that there are solutions, which are not consider herein 



LINEAR STABILITY FOR PERIODIC WAVES OF BOUSSINESQ AND KGZ 3 

the form y3(a; — ct), |c| G (—1, 1), which obey the equation (?2:a;[c^¥'+'/3" — (/5+/((/9)] =0, < x < T. 
whence, there exists a, C, so that 

c^Lf + (f" — Lf + f{tf) = Cx + a, < X < T 

By the periodicity of if, we conclude that C = and thus, we have a family of waves satisfying 

(3) f" - (1 - c^)^ + fi^) = a, < X < T. 

We now construct solutions of Q in various cases of interest, most notably p = 2 and p = 3. This 
material is not new, but in order to introduce the particular parametrization that is convenient 
for us, we include a sketch of the construction for completeness. 

2.1.1. The case p = 2. We will consider only the symmetric case a = 0. For the nonlinearity, 
/(n) = we have (with lo = 1 — c^) 

(4) -wip+ ^ip'^ + if" = 0. 
Therefore, 

(5) ^'2 _ ^^2 ^1^3^^ 

Hence the periodic solutions are given by the periodic trajectories H{(p, (p') = 6 of the Hamil- 
tonian vector field dH = where 

H{x,y) = y'^ + ^ - wx^. 

The level set H{x,y) = b contains two periodic trajectories if if > 0, 6 G (— |t/;^,0) and a unique 
periodic trajectory if 6 > 0. Here we consider the cases where 6 < and ipc > 0. To express (pc 
through elliptic functions, we denote hy ipi > ipo > the positive solutions of gP^ — — b = 0. 
Then ifo ^ 'Pc ^ and one can rewrite ([S]) as 

(6) if'^ = - 99o)((/3l - ^c){h^c + (pf) + ^l- "iuj). 

Introducing a new variable s G (0, 1) via Lpc = 1^0 + {^1 ~ 'Pq)s'^ ■, we transform ([g]) into 

where a, k, k' are positive constants (A;^ + k'"^ = 1) given by 

2 2ipi + ipQ - 'iuj 2 V^i - ^0 

a = — , k — 



12 2tpi + tfQ — 

Therefore 

(7) Vc{x) = <^o + {^i - ^o)cv?{ax] k). 

By the above formulas, 

'Pi — 'Po = 12a^K^ 
pi = 4a2(l + K?) + w 
ipQ = 4a2(l - 2^2) + w 



(8) 
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The fundamental period of the cnoidal wave Lp^. in ([T]) is 

9 T = ^^= ^ >^ , TG — ,00 

a \/w \ \/w 



Here and below, K{k) and E{k) denote the elliptic integrals of the first and second kind in a 
Legendre form. Further, we will use the following relations 

E{k) - (1 - k')K[k) E{k)-K{k) 
^^^) = ' ^ ■ 

Lemma 1. For any w > and T € (;^,oo), there is a constant b = b{w) such that the periodic 
traveling solution ^ has period T. The function b{w) is dijjerentiable. 

Proof. It is easily seen that the period T is a strictly increasing function of k: 
d . 4/- 5 T,., k(2k2 - 1)K{k) + 2(1 - k2 + k'^)K'{k) 



3 



2{l- + k^)E{k) + {1- t^){K'^ -2)K{k) ^ 

4k(1-k2)^1_k2 + ^43 

Given w and b in their range, consider the functions (foiw, b), ifi{w, c), k{w, b) and T{w, b) given 
by the formulas ^ and ([s]). We obtain 

dT _dTdk _ 1 dT d{k'^) 
~db ~ 'dk'dh ~ 2k~dk db 
Further, using that ^^Pq — = ~ wy'iiWe have 

db (2(^1 -\- ifQ — 3w)'^ 



{(fl - 2w(po){fi - '^w(fi){2(pi + (po- 3w) 



2 ■ 



We see that ^'^^^'^^ 7^ 0, whence the implicit function theorem implies the result. □ 

2.1.2. The case p = 3. We consider the "symmetric" case a = only, with f{u) = u^. Denote 
w = 1 — £ (0,1). Multiplying by if' and integrating implies 

4 

(10) = b + w^'^ 

Hence the periodic solutions are given by the periodic trajectories H{ip, (p') = b of the Hamiltonian 
vector field dH = where 

4 2 

H{x,y) = y + — - w — . 

Then there are two possibilities: 

• {outer case): for any 6 > the orbit defined by H{(p,ip') = b is periodic and oscillates 
around the eight-shaped loop H{(p, ip') = through the saddle at the origin. 

• {left and right cases): for any b G (— ^^^,0) there are two periodic orbits defined by 
H{ip, ip') = b (the left and right ones). These are located inside the eight-shaped loop and 
oscillate around the centers at (=F\A'^, 0), respectively. 
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We will consider the left and right cases of Duffing oscillator only. In these cases, denote by 
ipi > ifo > the positive roots of the quartic equation ^ — wz —6 = 0. Then, up to a 
translation, we obtain the respective explicit formulas 

/ N / N , / , N ,9 V^'i — V^'r\ 2o9? — 2w ipi ^ 2K(k) 
11) ^(z) = ^^idn{az;k), = ^ ° = , a = ^, T= 

Note that the fundamental period T may be written as follows 

/ X nn 2K(k) 2K(k)V2^^ rr. ( V^TT \ 

a Jw \ Jw / 



Note that this is a two parameter family of solutions, parametrized explicitly in this case by 
ipQ,ipi, although we shall need a different parametrization. In fact, we would like to think of 
this family as being parametrized (implicitly) in terms of T and w, where these two will be 
independent of each other. We have the following 

Lemma 2. For T > there is a constant b = b{w) such that the periodic traveling-wave 



solution (11) determined by H{ip,ip') = b{w) has a period T. In addition, the function b{w) is 



differentiable. 

For the proof, see Lemma 3.1 in 

2.2. Construction of the traveling wave solutions for the KGZ system. We are looking 
for T periodic traveling solutions of the Klein-Gordon-Zakharov system, ([2]). Thus, we take the 
ansatz u{t,x) = (pc{x — ct),n{t,x) = ipc{x — ct), where we take the speed c G (—1, 1). Plugging 
into Q, we obtain the following relation between ^/^ and (p 



Two integrations in x imply 



(c^ - l)ip{x) = '^f'^ix) + bx + a, 



for some constants a,b. By the periodicity we have that 6 = 0, whence V'(^) = — 2(i-c^" • 
simplicity, we shall consider the case a = only. That is 

Returning back to the other equation in Q and using the relation (13), we obtain the following 
equation for ipc 

(14) _(i_c2)^;' + <^^__^ = o<x<r. 
Thus, as in (1.1.2) after multiplying by ipc and integrating we get 

(15) ^? = ^ + ^ - ^ = 4^(4-^^ + ^w^l - ^t). 



w 



Denote by > 920 > the positive roots of the polynomial P{z) = z — Awz — Aw b. Then ( 15 ) 
can be written in the form 

(16) v'c=^{^l-'Pl){^l-^l)- 
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and the solution of (15) is given by 



(17) ipc{x) = ipidn{ax, k), 

where 



(18) ^l + V>l = ^^, « = ^, .^ = ^H^. 

2w iff 

Moreover, 

(19) {2-K^)ipl = 4w, a= I - AwH = Awy^l - iff. 

Since dn has fundamental period 2K{k), then the solution ipc has fundamental period T = 
In terms of k,w, given by 

(20) T = 2K{n)y/2 - k^^^, T e I = iV2^T^/I^^,oo) 



Lemma 3. For any w > and T > \f2TX\fw, there is a constant b = b{w) such that the periodic 
traveling solution JTa) has period T. 



Proof. The period T is a strictly increasing function of k: 



From (18) and (19), we have 

dT dTdK 1 dTdK? 



db dn db 2k dn db 
dK? dK? d(fi 16w^ 



db dipj db ipf{4w-2ipl)' 
The implicit function theorem then implies the result. 

□ 

3. Results 

3.1. Setting the linear stability problem for the Boussinesq equation. We now set up 

the linear stability/instability problem for ([T]). Set the ansatz u = (pc{x + ct) + v{t,x + ct) and 
ignore all terms O(u^). We get vu + 2cvtx + Mv = 0, where 

Mv = dtv - (1 - ^)dlv + U'i'fMxx 

Note that this operator M is not self-adjoint. However, if we introduce the variable z : Zx = v, 
we get the following linearized equation in terms of z, 

(21) zttx + 2cztxx + M{zx\ = 0. 
Note that M[zx] = dx[H[z]], where 

(22) H,z = d^z - (1 - c^)dlz + {f'{ipc)zx)x 

Thus, the linearized equation becomes dx[ztt + 2cztx + Hz] = 0. In our considerations, we say that 
the wave ipc is spectrally unstable, exactly when there is an exponentially growing mode, that is a 
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pair A G C : KA > and a T periodic function ^ G H^eA^^ T), so that d^iX^tp + 2c\^' + Hip] = 0. 
This of course imphes upon integration that for some constant a 

A^ + 2c\ip' + H^ = a. 

Integrating in [0,T] and taking into account that both ip'^H'ip are exact derivatives, impies that 

A*^ f"^ ih(x)dx ~ pjT ~ 

a = — '^^^ . Thus letting ijj := tp — implies that Jq ^{x)dx = and 

X^i) + 2cXtp + Hi> = 
These arguments motivate the following 

Definition 1. We say that the traveling wave ipc is spectrally /linearly unstable, if there exists an 
T periodic, mean value zero function G D{Hc) and X : KA > 0, so that 

(23) AV + 2cAV'' + Help = 0. 
The question for linear stability of equations in the form 

(24) Ztt + 2fizt^ + 'Hz = 0, 

or what is equivalent (at least in this case) to the solvability of 

(25) AV + 2/iAV'' + -HV' = 0. 

in Lq(0,T) = {/ G Lpg^(0, r) : f{x)dx = 0} has been addressed in a recent paper by the 
second and third author, [26]. Note that the self-adjoint operator T-L that appears in (22) is in 
the form 

n = -a,xa,, c = -dl + (i - c^) - f\^,). 

Here L is the ubiquitous standard second order Schrodinger operator, which appears in the lin- 
earization of the generalized KdV equation around its traveling wave solution ipc- This observation 
will be crucial for the spectral properties of the operator H as the properties of L are generally 
well-known, at least for the cases into consideration, p = 2, 3. 

3.2. Setting the linear stability problem for the KGZ system. We linearize the KGZ 
system as follows. We take u(t, x) = (^c(a^ — ct) + v{t^ x — ct) and n(t, x) = ipdx — ct) + h{t, x — ct) 
and ignore the contributions of all quadratic and higher order terms. We obtain the following 
linear system for the corrections v,h, 

vtt - 2cvtx - (1 - c'^)Vxx + V + IpcV + iPch = 
htt - 2chtx - (1 - c^)hxx - {Vcv)xx = 0. 

Further, introducing a new mean- value zero function z, so that h = Zx and w = 1 — (? . The 
second equation in (26) becomes 

dx[ztt - 2cZtx - (1 - C^)Zxx - i<^cZx)x] = 0, 

whence integrating in x yields 

Ztt - 2cZtx - (1 - c'^)Zxx - i'^cZx)x = a{t) 

for some function o(t). Observe however, that in our choice of z, we have required that 
J^^ z{t,x)dx = 0. Thus, integrating the last equation in [0, T] yields that all integrals on the left 
are zerc|^a(t) = 0, whence ztt — 2cztx — (1 — c'^)zxx — {^cZx)x = 0. We have shown that one can 
rewrite the linear stability problem (26) as follows 

(27) ^tt - 2c$tx + n$ = 0, 



(26) 



■^each term is either an exact derivative or Zu, which is mean value zero 
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where $ = ( ^ ) and 



(28) ^ - ( 5 H2 

where 

2 

H, = -(1 - c^)dl + 1 + V'c = -wdl + 1 - 
H2 = -wdl; 

Az = (fcz^; A*z = -{ipcz)x- 

Clearly, the operator T-L is self-adjoint, when considered over the domain 

Din) = H^[0,T] X H^[0,T] C L^[0,T] x L^[0,r] 

so that T-L : D{'H) — t- L^[0,T] x Lq[0,T]. Note that , when considering this spectral problem, our 
basic Hilbertian space will be L'^[0,T] x L^[0,T], instead of the usual L^[0,T] x L^[0,T]. 

3.3. Precise formulation of the main results. Our main results are described in the next 
Theorems. Note that in all of them, the issue is linear stability for the said families of spatially 
periodic solutions, when the perturbation is taken to be periodic with the same period as the 
underlying solution. 

The issue for stability/instability, when perturbations are taken to be with period of the form 
nT, n integer is a more complicated one. Clearly, the instability results continue to apply in 
this case, but it may very well be that some, previously stable waves (in the context of the same 
period perturbations) become unstable, when perturbed by nT periodic functions. 

Theorem 1. Let the nonlinearity in ([T]) has the form f{u) = v? . Then, the two parameter family 
of dnoidal solutions, described in (llll), are linearly stable, if and only if 



\c\ > \ -r^^^,Ke (0,1) 

where 

_ [4E{^) - [(2 - ^^)E{k) - 2(1 - k^)K{k)] 

^^""^ {2-k'^){E^{k)-{1-k^)K{k)) ' " 

where E{k), K{k.) are elliptic integrals of the first and second kind in a Legendre form. 

We now give a different formulation of the main result. Let T > v^tt. Then, the waves 



described in (11) are one parameter family of waves, having a fundamental period T, which can 



be parametrized by c : |c| < a/1 — (note that c, k are in one-to-one relation given by (12)) 



Now, Theorem [T] asserts that the stable waves in this family are exactly those with |c| > ct, 



where ct € (0, y 1 — w) is determined as follows. Let kt be the unique solution of 



K{k)V2 - k;V4 + M(k) = T 



Then, ct 



M(kt) 



4:+M(kt) ■ 

Our next theorem concerns the quadratic case 
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0.4 0.6 0.8 1.0 




, together with its terminal 



Figure 2. The blue line is the graph of the increasing function 
K{k)V2 - K^^/4: + M{k) : (0,1) The red line is V27r. Note that 

the range of this function is (v^tt, oo). 



Theorem 2. Let the nonlinearity in ([T]) has the form f[u) = \ 
are linearly stable, if and only if 



Then the periodic solutions 



c > 



Fin) 



4 + F(k) 



,K G (0,1) 



where 



2F(k) 



16Vl-«2+K4ii'2(K) 
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and 
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F{k) = 16K{k)[3E{k) + {k^ -2 + y^l - k"^ + k^)K{k)], 



1 



128 



dn 



0.5 ■ 



03 ^ 



07 



0.1 ^ 




0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3. Graph of the positive function 
value 2 



4+F(k) 



, together with its terminal 



An alternative formulation is the following. Fix T > 27r and consider the family of cnoidal 
solutions, described in ([T]). This is a one-parameter family of solution^ having fundamental 

period T and indexed by say c, where c : c < y^l — Then, the linearly stable waves with 
period T are exactly those, for which 



47r2 , , 

where ct £ (0, y^l — is determined as follows. Take kt to be the unique solution of the 
algebraic equation 



2K{k) \/1-k'^ + k^J 4 + F{k) = T. 



Then ct = , , . 

^ Y 4+F{kt) 

Remark: Using the results of Theorem [T] and Theorem [2j one can reconstruct the results on 
linear stability of the whole line waves, [26] 



(29) 



p + 1 



(1 - c^) 



1 

p-1 



sechp- 



Recall that the results of [26] predict linear stability if and only if |c| > 



^where n and c are related by ^ 
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Figure 4. The blue 



line is 

k^ + k\U + F{k) : (0,1) 



that the range of this function is (27r,oo). 



the 



graph of the function 
The red line is 2ti. Note 



and limK->i-(l 



Take p = 3. Then, the periodic waves described in (11) in the limit k 
whole line waves described in (29). Observe that since lim^-s-i- E{k) 



0, we conclude easily lim^^i- M{k) 
M{k) \/2 



1— correspond to the 
l,limK^i_ K{k) = oo 
4, whence 



lim 

K^l- 



Similarly, one can check that for p 



4 + M{k) 
2 (more precisely if f{u) 



\)), we have 



lim 

K^l- 



F(«) 



4 + F(k) 2 

Thus, we obtain the results in [26] for p = 2, 3 as a corollary. 
Our next result concerns the KGZ system 

Theorem 3. The KGZ system ([2) has a two-parameter family of traveling wave solutions 
{iPcipc); described in (13) and (17). These waves are stable, if and only if k, G {kq,1) and 

1 



1 > c > 



Vl + 4iV(K)' 



where kq = 0.937095... and the function N is defined in (73). 
If we take a limit as k — t- 1, we have 

1 V2 

iim 



v^l + 4iV(K) 2 ' 

Since k — )• 1 corresponds to the case T = oo or the case of the whole line, this allows us to 
conclude that the corresponding whole line solitons are stable, provided |c| > This was 

indeed the conclusion in [26], so we are able to deduce this result, as a consequence of Theorem 
|3J 

Once again, we provide an alternative interpretation of Theorem [3} Let T > be a fixed 
period. Then, there exists a one parameter family of periodic waves with period T, described in 



12 



SEVDZHAN HAKKAEV, MILENA STANISLAVOVA, AND ATANAS STEFANOV 
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0.9 ^ 



0.8 ^ 



0.7 



0.6 ^ 




Figure 5. The red line is the graph of the function 



green hne is the terminal value of 



V2 



^1+4:N{k) 



in (kq, 1). The 



(13), (17). More precisely, this family may be parameterizec^ by c with the following restrictions 

on c: if T < y^vr, then 1 > |c| > y^'l ~ ^j^i otherwise if T > \/27r, c G (—1,1). Then, the stable 
waves in this family are given by 

, , 1 

c > Ct 



V^l + AN{kt) ' 

where kt G («^0) 1) is found as the unique solution (see graph below) of the algebraic equation 



Vl + 4iV(K) 



T. 




Figure 6. Graph of the increasing function V-^(^) . i ) 

Note that its range is (0, oo). 



'The parameters k and c are related by ( 20 1 
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4. Preliminaries 

4.1. Linear stability theory for second order equations. In this section, we give a precise 
statements of tiie results of [_26j, concerning tlie linear stability of (24) or what is equivalent the 
solvability of (25). We assume the following about the self-adjoint operator %: 

(30) 



a{'H) = {-52} u {0} U a+{U),a+{U) C (a^, oo), a > 
Ucj) = -d^(j), dim[Ker{n + 6^)] = 1 
n^Jo = 0, dim[Ker(n)] = 1 
llV'oll = 1 



Next, we require that for all r >> 1 (note H + t > will be invertible) 

(31) {n + ry^'^d^in + t)-'I\ {% + ry^d^ e b{l^), 

Finally, we require 

(32) m = nfi. 

Note that the last identity ensures that H maps real-valued functions into real-valued functions. 
The following theorem, in a more general form, appears as Theorem 1 in |26j . 



Theorem 4. Let % he a self-adjoint operator on a Hilber space H. Assume that it satisfies the 
structural assumptions ( [SO] ), ([31]) as well as the the reality assumption (32). 

Then, i/('?^~^['(/'o], V'o) ^ 0, one has a solution to (25) for all values of fj, £ R^, that is instability 
in the sense of Definition^ Otherwise, supposing {'H~^['iP'q\,iI^q) < 0, 

• the problem (p5|) has solutions, if ijl satisfies the inequality 



(33) 



< |//| < 



= :/i*(^) 



• the problem (25) does not have solutions (i.e. stability), ifuj satisfies the reverse inequality 

(34) \ii\>iJ*{n) 

Remarks: 

(1) Theorem [4] appears in |26j as a result about the stability of (24), but we state it in its 
equivalent form for solvability of (25). 

(2) In the applications that we consider, we restrict our attention to the Hilbert space H = 
-Lq[0,T] or H = L'^[0,T] X Lq[0,T], depending on the situation that we are in. 



4.2. Spectral theory for the Schrodinger operators of Boussinesq waves. We review and 
state the main results regarding the spectral theory for the second order Schrodinger operators, 
arising in the linearization around Boussinesq waves. We consider again the cases p = 2 and 
p = 3 separately. 



4.2.1. The case p = 2. In this section, we present some spectral results for C that will be useful 
in the sequel. The first is a technical lemma, that will be used to establish the simplicity of the 
zero eigenvalue for 7i = —dxCdx. 



Lemma 4. (£"H,1) / 0. 
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Proof. In this case C = —d'^+w — ^pc and KerC = spamp'^. The spectral properties of the operator 
A = - 4(1 + k^) + 12Psn'^{y;k) in [0,2K{k)] are well-known [17j. The first three(simple) 
eigenvalues and corresponding eigenfunctions of A are 



^0 = - 2 - 2^/1 - k2 + 4k4 < 0, 



ipo{x) = dn{x; k)[1 - (1 + 2^2 - y/l - k"^ + 4:K'^)sn'^{x; k)] > 
m = 

V'i(x) = dn{x; K)sn{y; K)cn{ax; k) = \-^cn'^{y'-, 

= k2 - 2 + 2^/1 - k2 + 4^4 > 
^^2(3;) = dn(x; - (1 + 2^2 + ^/Y^^Kp^T^)sn^{y; k)]. 

Since the eigenvalues of C and A are related by = a'^fJ-n, it follows that the first three eigenvalues 
of the operator £, equipped with periodic boundary condition on [0, 2K{k)] are simple and 
Ao < 0, Ai = 0, A2 > 0. The corresponding eigenfunctions are ipo{ax),7pi{ax) = const.tp' and 
ip2{ax). Note that 1, (/9c -L KerC and 

(35) C{l)=w-^c 



and hence, we can take inverse in (35) 

(36) 1 = tt;£-H - £- Vc- 

Taking dot product with 1 yields the relation 

1) = -(1, !) + -(£- Vc,l) 
w w 

Differentiating ^ with respect to c, we get = 2apc: whence 

Entering this last formula in the expression for 1), 

(38) 1) = It + f 9c / V'cda; 

w 2cw V Jo 

Using (8) and that /q^^^"^ cn2(y; K)(iy = ^[E{k) - (1 - K2)/i:(K;)], we get 

(39) / (pcdx = 8a[3E{n) + (k^ - 2 + V^l - + k^)K{k)] = -F{k), 
Jo T 

where 

F{k) = IQK{k)[3E{h) + {k^ -2 + Vl - k2 + h^)K{k)]. 
Now, we need to compute 

dcFU) = F'U)- = -2cF'U) — . 

^ ' ^ ' dw dc ^ ' dw 

Thus, to compute we differentiate with respect to w the following relation 

(40) w^ = 16a^(l -K^ + K^) = 256 ^ — 
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obtained from (Isl) and (l9|). We get 
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(41) 
where 



dK 
dw 



G{k) 



128 



wT^G{k), 



1 

d[K^(K)(l-K'^+K'i)] 

dn 



From the above relations, ([9|) and (40), and we have 

(42) (£-il, 1) = - - T'^F'{k)G{k) = -[l- 16\/l - + k^K'^{k)F'{k)G{k)]. 

w w 

The expression in the brackets above is strictly positive (see the graphic below) which proves the 
lemma. □ 



035 ^ 



030 ^ 



0^5 ^ 



0.15 ^ 



Figure 7. Graph of the positive function [1 - 16\/l - k? + k'^K^{k)F'{k)G{k) 



4.2.2. The case p = 3. Consider 

(43) C = -dl + w-Sip"^. 



We use ((11)) to rewrite the operator C in an appropriate form. From the expression for (p{x) 
from (11) and the relations between the elliptic functions sn{x), cn{x) and dn{x), we obtain 

C = a^[-dl + Qk'^sn'^iy) -A-k^] 

where y = ax. 

It is well-known [17] that the first five eigenvalues of A = —dy + Qk"^ sn? {y , k) , with periodic 
boundary conditions on ^,AK{k)], where K{k) is the complete elliptic integral of the first kind, 
are simple. These eigenvalues, with their corresponding eigenfunctions are as follows 





= 2 + 2k^ - 2V1 


- fc2 + k^, 


^o(y) 


= i-{i + e-Vi- 


- k^ + k^)sn'^{y,k) 




= l + fc2, 




^i(y) 


= cn{y, k)dn{y, k) = 


sn'{y,k), 


V2 


= 1 + 4A;2, 




V'2(y) 


= sn(y, k)dn{y, k) = 


-cn'{y,k), 




= 4 + A;2, 




V'3(y) 


= sn{y, k)cn{y, k) = 


— k~'^dn'{y, k), 




= 2 + 2A;2 + 2VI 


- P + A:^, 


i'iiy) 


= 1 - (1 + P + - 


- k^ + k*)sv?{y,k) 
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It follows that the first three eigenvalues of the operator £, equipped with periodic boundary 
condition on [0, 2K{k)] (that is, in the case of left and right family), are simple and Aq = 0^(1^0 — 
1^3) < 0, Ai = a^(f3 — fa) = 0, A2 = a'^{i^4 — z^s) > 0. The corresponding eigenfunctions are 
Xo = ^0(03;), xi = f'{x),X2 = iJ4{ax). Thus, we have proved the following 



Proposition 1. The linear operator C defined by (43) has the following spectral properties: 

(i) The first three eigenvalues of C are simple. 

(ii) The second eigenvalue of C is Xi = 0, which is simple. 

(iii) The rest of the spectrum consists of a discrete set of eigenvalues, which are strictly positive. 
Next, we verify the following technical result. 

Lemma 5. The operator C verifies {C^^l, 1) 7^ 0. 



Proof. This statement was needed and proved in the work of Deconinck-Kapitula, [T^, but we 
repeat the short argument for completeness. First we will prove that for i ^ 0,4, {ipi, 1) = 0. 
Using the expression for A and ■04, we get 

12k; 

Uiiipi, 1) = 6K'^{sn^{y; K),iJi) = (1 - ip4,i^i). 

1/4 

It follows that for i / 4, 

(44) = (V'4,V'i) = 

Observe however that I'oi'i = 12k^, which means that (l — 7^ 0' whenever i ^ (since 

Ui / fo)- By (|44]), this implies that (V'i, 1) = 0,z / 0,4. 

From (Theorem 2.15, in |22j), the eigenfunctions of the operator C form an orthonormal basis 
of L^[0,T] and hence we compute {C^^l, 1) by expanding 1 in the eigenfunction expansion. Note 
that all terms corresponding to mean value zero eigenfunctions disappear (since {ipi, 1) = 0,i ^ 
0,4). Hence, the expansion for {C~^l, 1) has only two non-zero terms. More precisely, we have 



a^(z^o - 1^3)1 IxolP a^(z^4 - 1^3)11X2! 



(45) 



2 - 2VT )B3{k) (k2 - 2 + 2VI - k2 + 

where we have used the formulas 

2K{k) 2 

sn\y)dy = ^[K{k) - E{k)]; 

r2K{K) 9 



sn\y)dy = ^ [(2 + k^)K{k) - 2(1 + k')E{k) 
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and 



- k2 + ^4 _ 1 1 + ^2 _ ^1 _ ^2 + ^4 
5 A (k) H 5 



I -2 + E{k) 

- ^ % '-[K{k) - E{k)1 



(1 + k2 _ Vl - k2 + ^4)5 



[(2 + K2)i^(K) -2(1 + k2)^(k)]; 



Bi{K) = K{k)-^^ ^-^^ — '-[K{k) - E{k)] 



(1 + ^2 ^ _ ^2 ^ ^4)2 

3k4 



(2 + k^)K{k) - 2(1 + k;2)£;(k) 



+ 



From the graph of the function — , 
that (£^^1, 1) > and hence Lemma [s] is estabhshed. 



below, we reahze 
□ 




Figure 8. Graph of the positive function 



+ 



_B2M_ 



(k2-2-2v'1-k2+k'*)B3(k) {k2_2+2v^1-k2+k4)B4(k) 
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Lemma |4] and Lemma [5] allow us to verify the important property about the simplicity of the 
zero eigenvalue for the operator %c- 

Corollary 1. Let p = 2 or p = 2>. Then, the operator Tic = —dxCdx (corresponding to f{z) = 



,2 



z , ^) defined in (22) has zero as a simple eigenvalue in Lq[0,T], with an eigenfunction 



i'o = ^c - T k Re- 
proof. First (fc — )p f ipc{x)dx is easily seen to be an eigenfunction, since 

Ticifc - ^ I fc{x)dx] = -dxC[(f'c] = 0, 



since ip'^ is an eigenfunction for C. 

Regarding uniqueness, let / G Lg[0,T], so that Hcf = 0. It follows that 

^[/'] = c = const. 

Since Ker{C) = span{ip'^} and 1 _L ip'^, we can resolve the last equation as /' = cC~^l. Thus 

= (l,/) = c(l,£-il), 

whence c = 0, since 7^ by Lemma [sj It follows that /' is an eigenvector for C. Thus, 

/' = fi(p'^, by Propositionhl But then, / = id[ipc — f 'Pc{x)dx], since we are in the space Lq(0, T) 
and the uniqueness is established. □ 

4.3. Spectral theory for the Schrodinger operator T-L of the KGZ system. First, as in 
the Boussinesq case, we show that the operator T-L has a simple eigenvalue at zero. In addition, 
we identify the unique (up to a multiplicative constant) eigenfunction of H. Recall that in our 
considerations, we work with the space Lg[0,T], that is the second component will contain only 
functions with mean value zero. This proposition mirrors closely the corresponding statement of 
Proposition 8 in |26j . with a few notable differences. 



Proposition 2. The self-adjoint operator T-L introduced in (28) has an eigenvalue at zero, which 



is simple. In addition, the unique (up to a multiplicative constant) eigenfunction is given by 



\ 2w 



Proof. Let ^ ^ ^ be an eigenvector corresponding to a zero eigenvalue, that is ?^ ^ ^ ) = 0. In 
other words, 

(46) I -^f" + f-^f + ^9' = 

Integrating the second equation in x implies that for some constant cq, we have 

g = h Co, 

w 



whence the equation for / becomes 

(47) _y,/" + /_|^/ + co<^ = 

We will show that cq = and then / = d(p'^ for some constant d. To that end, recall the defining 



equation for tpc, namely (14) and differentiate it with respect to x. We get 



(48) _^^;f + ^;_g^; = o. 



LINEAR STABILITY FOR PERIODIC WAVES OF BOUSSINESQ AND KGZ 19 

Following the usual analogy with the KdV equation, we introduce the second order differential 
operator 



" ^ 2w 

Using that dn'^ + K^sn? = 1 and wa^ = gz^j get 

'2 



C = —wd^ + 1 — "^{1 — K^sn'^iax, k)) 
= —wd^ + 1 — 6wa^ + G'wa'^K?s'n?{ax, k) 



wa^{—dy + 6sn^(y, k) — — 4), 



where y = ax. It follows that the first three eigenvalues of the operator £, equipped with periodic 
boundary condition on [0,2K(A;)], are simple and Aq = wa^{uQ — v^) < 0, Ai = wa'^{v^ — t's) = 
0, A2 = wa'^{vi — 1/3) > 0. The corresponding eigenfunctions are xo = V'o(a2;),xi = ^'{x),X2 = 
V'4(a;x), where Vi and tpi are given in (3.2.2). 

In particular, the kernel of C is spanned by (p'^, i.e. Ker{C) = span[ip'^]. Going back to (47), 
we can rewrite it as 

C[f] + CQifc = 0. 

Note that all solutions to this last equation are given by 

f = dip' - coC^lipc], 

where d is an arbitrary scalar, since ipc -L span[ip'^] = Ker[C). Plugging this last formula in the 
equation for g yields 

g = h Co = [d^p -cqL 99 ) + Co = (f^p + cq h 1 

w w w \ w J 

Integrating the last expression in [0, T] and using the periodicity yields 
(49) co(^^^^^^+t)=0. 



Thus, if we verify that {(p,C ^[(p]) 7^ —Tw, we would conclude from (49) that cq = 0, whence 
/ = d(f'. Furthermore, g' = —^ipp', whence 

g = —d h const. 

2w 

Recall however that the const, in the formula above is uniquely determined by the fact that g 
has mean value zero (i.e. g £ Lq[0,T]), whence 



9 = d[ + 



2w 2Tw I 

Thus, Proposition [2] is established modulo the following 

Fact: {(pc C~^[(pc]) > so in particular {(p,C^^[(p]) > —Tw 

From (14), we have ip^ = —2'w'^ip'^ + 2wipc and thus 

3 

(50) Lip = -w<p'^ + ipc- —^l = 2wip'^ - 2lpc- 
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o3 



On the other hand, differentiating —w\p^ + wip, 
w) and using (50) to express results in 



0, with respect to w (and dividing by 



C 



difc 
dw 



//I 1 1 ^ 

2(/7 ipc = —^c H >C93c, 

www 



Taking £ ^ in the last identity yields 



'Oc = w-^ - ipc- 



Since /q^^"^ dn^d/, K)dy = E{k), we compute ip^dx = ^^^E{h)K{k). Therefore 



dw 



w^ A6w^ 
—o,„\—— 



E{k)K{k)] 



IGw" 



E{k)K{k) 



I (IK, aw 



To compute J^, note that from (18), we have 2wa = ipi, whence 



(51) 



'iw{2 - k^)K'^{h) = T^. 



Differentiating (51) respect to we get 
(52) 



dK 
dw 



Thus, using (51 ), 



8w^ d 
T dK 



-wT 



\K(n)E(K)] (^-^^)^^(^) 
^ ' ^ ' >'1[{2-k^)K^{k)] 

1[K{k)E{k)] 



l[{2 - K^)mK)] ^ 

Looking at the graph below, we realize that since | = lim^-^o 2-^^^^'^'^^^'^^^ 
we have that 



:[(2-k^)K^(k)] 



> 2 



f^[K{K)E(K)] 



{£. VcV^c) > -\wT, 



which establishes the claim. 



□ 



The next thing one needs to establish, in order to apply Theorem [4] is that the operator Ti 
for the KGZ system ( defined in (28) ) has a simple negative eigenvalue. This result should be 
compared with the corresponding statement in Proposition 9 in [26j for the whole line case. 



Proposition 3. The operator Ti, defined in (28) has a simple negative eigenvalue. 



Proof. Consider the eigenvalue problem in the form 



n 



f 
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0.338 h 



0J36 - 



0334 - 



0.332 - 



0.330 - 



0.328 - 




Figure 9. Graph of the function 2- 



f,mK)E(^)] 



for some o G (0, oo). As in Proposition [2| this can be rewritten a^ 
(53) 



Form the second equation, we may resolve for g 



(54) 



g = {a'-wdl)-^dM]- 



This last formula requires a bit of justification, but basically since dx [-pf] is guranteed to have 
mean value zero, it suffices to define (a^ — wd'^)~'^ (where > 0, w > 0) on L2[0,T] by 



(a^ 



a„.e 



E 



+ 47r2 



_g27rm|^ 



W 



rp2 



whence the formula for g in (54) makes sense. In fact, LQ[0,r] is invariant under the action of 
and (a^ — wd'^)^^ and hence g E L^[0, T]. In fact, we use ([54]) to deduce the following formula 
for g' 



9 



w w 



We plug this in the first equation of (53) to obtain the following equation for / 
(55) - wf" + (1 + a^)f - ^/ + ^[^c(a2 - wdly'^J] = 0. 



Introduce a one-parameter family of self-adjoint operators 



Ma := -wdl + (1 + 



+ 



In order to finish the proof of Proposition |3j we need to establish that there is unique oq > 0; 
so that the operator Mag has an eigenvalue zero and such an eigenvalue is simple. We shall first 
show that there exists ao > 0, so that Mag has an eigenvalue at zero and then we show that this 
eigenvalue zero is simple for Mag. 



'^Note that the second equation requires g{x)dx = 0. This means in particular that the statement for 
existence of negative eigenvalue is invalid, unless the second component of the Hilbert space is Lo[0,T] 
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To that end, we wiU estabhsh the foUowing 
Claim: For a > 6 > 0, we have > Mf, + (a^ _ p)id > Mh- 

Assuming the validity of the Claim, let us finish the proof of Proposition [sj Define A(a) to be 
the minimal eigenvalue for Ma, that is 

A(a) := inf{A : A G a(M„)} = Jnf^(M,/,/). 

Clearly, the function a — )• A(a) is continuous (in fact, more generally, o — )■ Ma is continuous as a 
function from — )• B(L'^[0,T]))). Moreover, as a consequence of the Claim, A(a) is a strictly 
increasing function of its argument. In order to show the existence of oq, it will suffice to show 
that a — )■ A(a) changes sign in [0, oo) (and hence vanishes at some ao > 0). But at a = 0, we have 



Mo = -wdl + l-^ = C, 



2w 

that was considered before. Since we have checked that C has a (simple) negative eigenvalue, say 
—(5^ it follows that A(0) < 0. On the other hand, by the claim 

Ma>Mo + a^Id = C + a^Id. 

In particular, for every a > 6, we have Ma > (a^ — 6'^)Id, whence A(a) > — 5^ > 0. Thus, for 
any b > 6, the function A(a) changes sign (exactly once) in the interval (0,6). We have shown 
that there is oq : A(ao) = 0. 

Now, we have to show the second part of the proposition, namely that is an isolated eigenvalue 
for Maa- Let cpo be the eigenvector for the simple negative eigenvalue for Mq = C. Note that 
by Proposition jij the second eigenvalue of C is zero, which means that Mo\^^± = >C|^^± > 0. In 
addition, by the Claim, we have that Ma^ > Mq + a^Id = C + a^Id, and thus, by the Courant 
maxmin principle for the second eigenvalue 

Al(Maf,) = sup inf (Ma^UjU) > flg + ™f {Cu,u) > Qq > 0. 
z^Q u±z:\\u\\ = l M±</)o:||u|| = l 

It follows that Ao(Ma(,) = 0, while Ai(Mao) > 0, which shows the simplicity of the zero eigenvalue 
for Mag. It now remains to establish the Claim. 

Proof. By the form of the operators Ma, it suffices to show for all trig polynomials / S L^[0,T] 
that for a > b, 

a'{^,{a^-wdl)-'[^J],f)>b'{^,{b'-wdl)-'[^J],f) 

But letting an be the Fourier coefficients of (pcf, that is ipcf = X/n Qn ^ or 
«n = /q^ (/Jc(2;)/(x)e~^™T(i2;, we see that 



b'{Mb'-wdl)-'[^J],f) = y^———^\an\'< 

IT- 
,2 



where we have used that j- < 5-, whenever w > a > b > {). □ 

□ 
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5. Linear stability for Boussinesq: proof of Theorem [T] and Theorem [2] 

Now that we have the sohitions, we need to check that the operator He satisfies the requirements 
in Theorem |4| after which, we need to compute the index ^*{Hc). We collect the needed results 
in the following propositions. 



Proposition 4. (Spectral properties of He) The operator He introduced in (22), satisfies (30), 
([32]) /or p = 2, 3. 

Our next result gives a precise formula for the index ^*{Hc)- 



Proposition 5. We have 
l)forp = 2 



F{k), 



where F{k) is defined in (61) below. 
2) for p = 3. 



lJ*{Hc 



/w 
2\ 



K(k) 



[{2- k'^)E{k)-2{1- k'^)K{k)] 



{2 - h'^){E'^{k) - {I - k'^)K{k)) 



Note: The function under the square root is positive for all values of k : < k < 1. 
We now finish the proof of Theorem [T| based on the results of Propositions |4] and [5} 
Let p = 3. We apply Theorem[4| whence we get stability, provided |c| > fi*{H). Thus, we need 
to resolve the inequality 



c > 



where we have taken 



M{k) :-- 



K{k) 



[{2- k^)E{k)-2{1- t^)K{K) 



{2-k?)[E'^{k)-{1-k'^)K[k)) 
we obtain for the interval of the stable speeds 



c > 



M(k) 



4 + M(k) 



as stated in Theorem [TJ 

Similarly, for p = 2, we have according to Proposition [5] that 



\c\ > ^ \/F{k) 



whence we conclude similarly that 



c > 



F(k) 



, 4 + F(k) 

is a necessary and sufficient condition for stability of the corresponding traveling wave. 
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5.1. Proof of Proposition |4| We note that standard arguments imply the vahdity of (31 ), since 
Ti is fourth order operator. The reality condition (32 ) is also trivially satisfied as all our potentials 
are real- valued. The hard to check condition is (30). We need to check that the operator Ti has 
a simple eigenvalue at zero. This was indeed the conclusion of Corollary [TJ 

Thus, it remains to show that the operator T-L defined in ( 22 ) has a simple negative eigenvalue 
and Proposition|4]will follow. This is a non-trivial fact. Interestingly enough, this was needed (and 
proved in our paper |[19j ) , when we have considered the transverse instability of the same spatially 
periodic waves in the Kadomtsev-Petviashvili (KP) and modified KP models, that is exactly for 
the operators considered here, corresponding to the cases /(u) = \ and f{u) = u^. More 
precisely, we have derived a necessary condition, so that the first two eigenvalues oiTi = —dxCdx 
satisfy 

Ao(^) < Ai(?^) = 0. 

which was verified for both p = 2, 3. This is exactly what is needed here. The interested reader 
may consult Section 4.1 in jl9| for full and complete proof of these facts. 



5.2. Proof of Proposition [5} 

5.2.1. The case p = 2. We compute the index of stability fi*{'H). We have 



(56) 

where A 
Hence 

(57) 

Thus 

(58) 



1 



'T Jo ^cdx. Let f : T-L[f] = (p'^. It follows that —Cf = (pc + b, for some constant b. 



f = c-^^c + bc-h. 



;{C ^iPcV'c) + b{C 4, 



^^ + ^||2N-"^^/ ||^^ + ^||2 

From (|57|), we have that = -(/', 1) = {C~^Pc, 1) + b{C-^l, 1), whence 



(59) 



b 



Combining (37), (56), and (59) yields 



From ^ after integrating 



1 



\v^c + A\\2 
1 



(£-11,1) 
1 d 



{(fc + AW^ yicdc 



(fj.dx = 2w ifcdx 



(£-Vc,i)(/:~H,(^e) 

(£-11,1) 

(£-Vc,l)(£-^l,V^c) 
(£-11,1) 



2w 



F{k), 



we use (41 ) to see 

1 d 



(60) 



4c dc 



T 



[-F{k) - 2mK\K)F'{K)G{K){l - k2 + K^)] . 



LINEAR STABILITY FOR PERIODIC WAVES OF BOUSSINESQ AND KGZ 



25 



Next, we use (37), (39), (41) and (40) to compute 



1 



2c 



(£-Vc, 1) = i^{dc / Mx) = -wT''F'{k)G{k) 



256 



[K\k)F'{k)G{k){1- + K^)] 

Finally, using the formulas for ip'^dx, ipcdx allows us to find 

2F{n) 



w 



\Wc + A\\'^ 
Putting all formulas together 



16V1-k2+k4k2(k) 



T 



W 



2F{k) 



-Fin) - 256A-(«)F»0(«)(1 - «^ + k') '"^'^'^'Mi^ - + 



1 - 16V1 -K^ + K^K^{K)F'{K)G{t 



2F{k) 



Thus, if we assign the function 
(61) 

F{k) + '. 

we get {7i~^ip' ji'') = — ^p^^-^ ■ Thus, the index formula holds as stated in Proposition jHj namely 



w 



5.2.2. The case p = 3. In this section, we compute the index of stability. For that we need to 
consider first 

(?^"Vo,V'o) = II,. , ^ii? (^~Vc^V^c). 



\^c + A\\^ 



where A = —ji ipcdx = — 2k(k) • Thus, we need to compute % ""^[i/^c]. Let f T-L[f] = <~p'c- It 
follows that —Cf = ipc + b, for some constant b. We conclude that 

-f = c-'ip, + bc-H. 

Note that C^^l is well-defined, since 1 _L c/?^, which spans Ker{C). Thus 
Differentiating ^ (with a = 0) with respect to c yields 



(62) (?^-Vo,V'o 



(63) 



From (11), we get (p'^dx = ^^('^^^in) 



(64) (^-Vc,^c) = ^(^,^c 



4c 



]_dApc 

2c dc 



ifldx] 



4 E'^{k) - (1 - K2)i^2(^) 



T 



k(1 — K^) 



26 SEVDZHAN HAKKAEV, MILENA STANISLAVOVA, AND ATANAS STEFANOV 

In addition, note that 



{C ^l,(pc) = {C Vc, l) = ^9c[/ (pc{x;K)dx] =0. 







and hence 



From the relations (11), we have 



w 



^ AK'^{k){2- k'^) 



2 ^1 r2 

which after differentiating with respect to w allows us to express 

dK r2 1 



Thus 



T{E^{k)-{1-k^)K^{k)) 1 



2 

Using that ^^^^ = ^'''"•'~(^l|!2)^^''^ , we obtain 



(65) 



{C ipc, ^c) 



1 E\k) - (1 - k')K\k) 



a {2 - k'^)E{h) -2{l- k?)K{k) a 
Since Jq^^'^^ dn{y; K)dy = vr, we get 



B{k). 



(66) 

From the above relations and ( |62[ ), we get 

(67) (H-Vo,V'o) 
Thus, 



K{k) 



aC{K) 



2-k' 
wC{k\ 



-B(k) 



1 



C{k) 



{2-k?)B{k) 



'W 



4E{k) 



K(k) 



[{2- k?)E{k)-2{1- k'^)K{k)] 



{2 - k^){E2{k) - (1 - k'^)K{k)) 
which is exactly the claim of Proposition [5] 

6. Linear stability of KGZ: proof of Theorem [3] 



We have already checked the conditions on the operator H, defined in (28) in Section 4.3 



Namely, we have established the simplicity of the eigenvalue at zero in Proposition [2] and then 
we have verified the existence and simplicity of a single negative eigenvalue. It now remains to 
compute the index ^*{%), after which, we obtain a characterization of the linear stability by 
Theorem |4j namely |c| > 



LINEAR STABILITY FOR PERIODIC WAVES OF BOUSSINESQ AND KGZ 



27 



Proposition 6. For k G (0, kq), kq = 0.937095 . . (?^" Vq, V'o) < 0. For k G (kq, 1), 
(T^-Vo.V'o) > 0, and 



where N is defined in (73). Note that, N{k) > 0, k E {kq, 1). 

Assuming the validity of Proposition [6| we now finish the proof of Theorem [3j To that end, 
observe that since {'H~^iI)q,iI)'q) < 0, k G (0, kq), we have instabihty whenever k G (0,fi;o)- For 
K G (kq, 1), we need to solve the inequality 



1 > c > 



w 



which results in the following necessary and sufficient condition for linear stability 

1 



1 > c > 



, K G (ko, 1). 



This is exactly the statement of Theorem [3j 

6.1. Proof of Proposition joj Now we will estimate the index of stability {'H~^tlj'Q,ip'Q), 

1 CT ,„2, 



ipo = m 



+ 



2w ^ 2wT 



where m is so that 



1. Thus, we need to compute T-l~ 



\2wl ■ 



We have 



Integrating once in the second equation yields 



)' 



(68) 



9 



w 



where ci is a constant of integration and needs to be determined. The first equation becomes 



2w 2w'^ w 



or 

(69) 



2u'2 



+ 



w 



It 



On the other hand, taking derivative with respect to w in (14) yields 
(70) 



From m'dv and (70), we have 



c 



dipc 



+ 



dw 2w^ 

d^c 
dw 



ci 

— ( 
w 
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and hence 



(71) 



difc ci 1 dip ci 

T = >C yjc = (l-ci)— H Lfc. 

aw w aw w 



Plugging this in ( 68 ) and integrating, we get 



1 f^,.2 



(72) 



ci 



1 rT,„2. 



T + 1 £- Vc> T + f,^ ^/^dx - 1 vldx ■ 



Using that j^^'^'^ ipl = ^E{k)K{k) and (pi), we have 



ipc-j—dx 
att; 



1 d 

2 ./o 



iiw 



2E{k)K{k) 



2dw' T 
{2-n^)K\K)i:[E{K)K{n)] 
l[{2 - n^)K\K)] 



The above formula allow us to express ci as a function of k only 



ci 



2^(K)j^(K)^[(2 - k^)K\k)] - 2(2 - ^'')K\k)1[E{k)K{k)] 

(2 - k2)K2(^) J.[(2 _ «2)^2(^)] _ 2(2 - K2)i^2(^) 

(2 - k2)^2(^) _ 8(;l _ ^2^e{k)K{k) + 2(1 - K^){2 - k'^)K^k) 
2(2 - k'^)^E{k)K{k) - 2(1 - k2)(2 - K2)i^2(^) _ 2(2 _ k^)E^{k) 



Now, 



' 2W' 



From (14) and the expression for / and g' , we get 



((^',',/) = ^(2J2- J3-J4 + J5), {g',^' 



1 



(Ji + J2-J3- Ja), 
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where 
h 

Ji 
J2 

J3 



^(''^ . 4 - 1 - A-2 

dn\y,K)dy = -——E{k) - -^K{k)- 

6 6 



2-^2 



Tci 



iVcVc) = 8ciE{k)K{k); 



T(l-ci)(y,fe,yg) 
2u;2 



8(1 -ci) 



2-^2 1 " j^[(2-Ac2)i^2(;^)](^ [2-^2 i_ 



J4 = 



J5 



^(V.„v^,)=32ci^. 



V^c, = 8(1 - ci) 

u; dw 



2K{k)E{k) 



{2-n'')K\K)i,[E{K)K{n)] 



dK 



[(2 - k'^)K^k)] 



In addition, we have 
1 



2«) 



l' = (<^c,¥''c) + |^(¥'c,<^c)-^ 



ifldx 



Using that dn'{y) = —K'^sn{y)cn{y) and sn'^{y) + cn'^{y) = 1, we get 



8w K{k) 



/ sn^{y,K)dy- / sn'^{y,K)dy 
Jo Jo 

[2(2 - k'^)E{k) - 4(1 - k^)K{k)]. 



,-2K{k) 



3T 2 - k2 



and 



w r 8i^(K) 

T [3(2-^2) 
16u; 



[2(2 - _ 4(1 _ k'^)K{k)] + 



2-k2 ^"^"3^ 



16 ^2 



T 



£;(K)^r(/«) 



1 - 



K\k) 



1 



r^'(«) 



2-/c2 ^ ^ 2-/^2" 
Combining the above relations yields 

-III ii\ Ji + 3J2 — 2 J3 — 2 J4 + J5 



(73) 



(H-Vo,V'o> = 



iV(Av) 



16u) 



E{n)K{K)-^^,K^{K)-^,E^{K) 



w 
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